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ELLIPTIC BOUNDARY VALUE PROBLEMS 
FOR THE INHOMOGENEOUS LAPLACE EQUATION ON 

BOUNDED DOMAINS 

STEVEN G. KRANTZ AND SONG-YING LI 



g 

Abstract. Elliptic estimates in Hardy classes are proved on domains with mini- 
mally smooth boundary. The methodology is different from the original methods 
^ ■ of Chang/Krantz/Stein. 



1. Introduction 



Let D be a bounded Lipschitz domain in W 1 . In [CKS], Chang, Krantz and Stein 
introduced certain real variable Hardy spaces. They defined two Hardy spaces h p (D) 
and h p (D), < p < 1. We say that a distribution / lies in h p (D) if it has an extension 
E(f) G h p (M. n ); and we say that a distribution g lies in h p z {D) if there is an extension 
E z (g) G h p (R n ) so that E z (f) =0onI n \D. Here h p (R n ) is a local version, due to 
Goldberg [G], of the classical Hardy spaces. 
, Since D is bounded, it is known from Miyachi [M] that the Hardy space h p (D) 

can be identified with the subspace of distributions / G T>(D) such that the radial 
maximal function f + (x) = fljj(x) G L P (D). Here 



™-,{|/ D «-.>/«* 



: < t < 8(x) 



where 5(x) = dis(x, D c ) with D c = R n \ D and is a fixed function such that 
G C^{B n ) (B n the unit ball in W 1 ), > 0, and f Rn <p(x)dx = 1. We set <j> t {x) = 



t- n (f)(x/t). 

^ ■ Let G(f) be the solution of the Dirichlet boundary problem: 



'1.1) Am = / in D, and u 



= 0; 

dD 



also let N(f) be the solution of the Neumann boundary value problem: 

(1.2) Au = f in D, and — =0, 

dv an 
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where v(x) denotes the outward unit normal to the boundary dD at x, and = V-za 
In [CKS], Chang, Krantz and Stein extended a classical theorem in IR n to a 
smoothly bounded domain in M. n and they proved the following theorem. 

THEOREM 1.1. Let D be a bounded domain in R n with C°° boundary. Let 
<p < 1. Then 



d 2 G{f) 



dxidxj 



+ 

h v r (D) 



d 2 N(f) 



dxidxj 



<C p (D)\\f\\wD). 

h p r {D) 



The question of what is the minimum smoothness condition on dD so that the 
above theorem remains true when D is a bounded domain in W 1 is still open. Cer- 
tainly one may ask: Does the smoothness of dD in Theorem 1.1 depends on pi In 
[CKS], it was conjectured that Theorem 1.1 remains true if D has C 1 ^ boundary. 

The Dirichlet problem for the inhomogeneous Laplace equation (1.1) has been 
studied by many authors (see, for example, [AND], [GT], [FKP], [LiM], and [Ken2]). 
Fix < a < oo and 1 < p < oo. When D is a smoothly bounded (C°°) domain in 
R™, Calderon-Zygmund theory shows (see [AND]) that if / G W*(D) (the Sobolev 
space), then there exists a unique u which solves (1.1) with 



(!-3) h\\wi +a (D) < C\\f\\ W p {D) , Kp<oo. 



In [Dah], B. Dahlberg constructed a bounded Lipschitz domain DcK 2 , and an 
/ G C°°(D) so that the regularity (1.3) fails when a — 0. In [JcK], D. Jerison and 
C. Kenig constructed a bounded domain in M 2 with C 1 boundary and a function 
/ G C°°(D) so that (1.3) fails for p = 1 with a = 0. Since C°°(D) C h\(D), we see 
that the aforementioned conjecture is not true when p — 1. 

/,From the definition of p-atom (see [CKS]), one can see that the necessary order 
of cancellation in an atom depends not only on p but also on n. Since the dual of H p 
is A n (i/ p _i) (the homogeneous Zygmund class), it seems that a reasonable necessary 
requirement on the smoothness of dD so that Theorem 1.1 remains true will be that 
the boundary is C n ^^ p ^ when p < n/(n + 1). 

The primary purpose of the present paper is to prove that this last suggested 
necessary condition is also sufficient. We will work on domains with at least C 2 
boundary. As we shall see, our results will be o interest for p small: < p < nj (n+2). 
We shall leave the case when D has C q boundary with 1 < q < 2 for a future paper. 

For < p < oo and any e > 0, we let 

(1.4) (x(p,e) = max{2 + e,n(l/p - 1) + e}. 



The first theorem we propose to prove is: 
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THEOREM 1.2. Let < p < oo and let D be a bounded domain in W 1 with C a{p ' e) 
boundary. Then 



d 2 G(f) 



dxidxj 



+ 



h p r (D) 



d 2 N{f) 



dxidxj 



h p r (D) 



<C p ^{D)\\f\\ K{D) 



for any e > 0. 



Note 1: From the proof of Theorem 1.2, it may be seen that we need such an e > 
only when n(l/p — 1) is an integer. 

We shall show by example that Theorem 1.2 fails for some domain D having only 

C n(i/ P -i)-e boundary. 

The proof of Theorem 1.1 given by Chang, Krantz and Stein is based on mapping 
D to the model domain R". Our approach for proving Theorem 1.2 will be based 
instead on the machinery connected with the fundamental properties of the Green's 
function for —A in D. The properties of the Green's function have, historically, 
played a crucial role in solving the Laplace equation. M. Griiter and K.-O. Widman 
[GW] as well as E. Fabes and W. Stroock [FaS] studied the Green's function in a 
Lipschitz domain in R n . They gave the basic estimates for the Green's function and 
its first derivative. In order to prove Theorem 1.2, we need the asymptotic behavior 
of G and its higher derivatives near the boundary. The secondary purpose of the 
present paper is to estimate the Green's function and its derivatives of all orders (see 
Theorem 2.2.) We believe that the properties of the Green's function that are derived 
in this paper will be useful in other contexts as well. 

Note 2: From our proof of Theorem 1.2, we have that 

d 2 G{f) 



dxidxj 



= r 2 (/) + c/ 2 (/), 



where r 2 is bounded from h p (D) to h p (D) for D with Lipschitz boundary. The 
operator U 2 : h p z (D) — * h p (D) is bounded, where h p (D) is the subspace of h p (D) 
consisting of all harmonic functions on D. 

^From our estimates on the higher derivatives of the Green's function, we are 
able to consider the elliptic boundary value problems (1.1) and (1.2) with / lying 
in the Hardy-Sobolev space. For this purpose, let us now introduce the definition of 
Hardy-Sobolev space. 

Let k be a non-negative integer and < p < oo. We let h k z ,p (D) denote the 
space of all measurable functions / with the weak derivative V / G h p (D) for all 
< £ < k; here V k f denotes all k th derivatives of /. We say that a measurable 
function / on D belongs to h k r > p {D) if V*/ G h p r (D) for all < £ < k. It is obvious 
that h k /(D) C h k > p (D). We shall prove the following: 
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THEOREM 1.3. Let < p < 1 and let k he a non-negative integer with k < 
n(l/p — 1). Let D be a bounded domain in R n with C z+[n{l/p - 1)] boundary. Then 



d 2 G{f) 



dxidxj 



+ 



h k r ' P {D) 



d 2 N(f) 



dxjdxj 



h k r ' v {D) 



< C Pik . n (D)\\f\\ h k, P{D) . 



The paper is organized as follows. In Section 2, we estimate the behavior of 
the higher derivatives of the Green's function near the boundary. The results in 
this section are the key to the rest of the paper. In Section 3, we use the atomic 
decomposition theorem in [CKS] to reduce the proof of Theorem 1.2 to the study 
of a single atom. In Section 4, we prove Theorems 1.2 and 1.3 for the Dirichlet 
boundary value problem. In Section 5, Theorems 1.2 and 1.3 for the Neumann 
boundary problem are proved and the 'Green's function' for the Neumann boundary 
value problem is studied. In Section 6, we shall construct examples that show that 
Theorem 1.2 is reasonably sharp. 

The authors wish to thank Jiaping Wang for a useful conversation on the Green's 
function. 



2. Estimate derivatives of Green function 
It is well-known that the fundamental solution for the Laplacian A in 



is 



r(,-») = r(|x-»|) = {^ o - | 5 1 ^ 



y\ 2 - n , ifn>2 
if n = 2 



where u n denotes the surface area of the unit sphere in W 1 (see [KR]). 

Let D be a bounded C 1 domain in IR n . For each y e D, we let U(-,y) be the 
solution of the Dirichlet problem: 

(2.1) AU(; y) = in D, U(x, y) = T(y -x), for x e 3D. 
Then the Green's function for A on D is 

(2.2) G(x,y) = T(x-y)-U(x,y). 
The Dirichlet problem (1.1) has a unique solution 

(2.3) u(x) = G(f)(x)= [ G(x,y)f(y)dy. 

Now let V(-,y) be the solution of the following Neumann problem: 



(2.4) 

We let 
(2.5) 



AV = Vo\{D)~ 1 in D, 



d 



V 



d 



dv dv 
N(x,y) =r(x-y)-V(x,y) 



F(y - •) on 3D. 
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It is easy to show that the Neumann problem (1.2) has a unique solution, up to an 
additive constant, given by the formula 

(2.6) u{x) = N(f)(x) = I N(x,y)f(y)dy. 

JD 

The main purpose of this section is to study the basic properties of U(x,y). We 
shall derive information about the asymptotic behavior of U (x, y) when x, y are near 
the boundary dD. For convenience, we will always assume that n > 2 (the case n = 2 
is similar, but the details of the formula are different.) Similar results for V(x, y) will 
be obtained in Section 5. 

We will need the following proposition that is due to M. Griiter and K.-O. Wid- 
man [GW]. 

Proposition 2.1. Let D be a bounded domain in M. n satisfying the uniform exterior 
sphere condition (i.e., each boundary point of D has an exterior osculating sphere of 
uniform size). Then the Poisson kernel 

satisfies the estimate: 

(2-7) 0<P(x,y)<^- 

\x — y\ n 

for all x G D and y G dD. [Note that, on a C 2 domain, these last two expressions 
are known to be comparable — see [KR]]. 

The main result of this section is to prove the following theorem about U. 

THEOREM 2.2. Let k > 2 be an integer. Let D be a bounded domain with C k °' e 
boundary. Then for any multi-indices a and (5 with \a\ = k < k and \(3\ — £ < k , 
we have 



(2.8) 



d k+£ U(x,y) 



(\x-y\ + 5(y)) n+k +^ 2 



dx@dy a 

for all x, y G D and any e > 0. 

Proof. If we write /(*) = -t 1 -"/ 2 for all t > 0, then T(x -y) = f(\x - y\ 2 ). Now 

where P(x, z) is the Poisson kernel satisfying (2.7) for all x G D and z G dD (since 
D has at least C 2 boundary) Thus 

P.10) ^=*gf»>, xedD , yeD . 
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For any fixed xq, y G D , we let R = \xq — y j/2. If R < 166 (y), then (2.8) holds 
with (x, y) = (xq, y) by the maximum principle and Proposition 2.1. We now assume 
that R > 166(y), and we let eo = 6(x ) n+£+k 6( y y) n+2k . For any multi-index a with 
\a\ = k, we consider the function 

d k f(\x-z\ 2 + e ) 



f a , y {\x-y\ 2 + e ) = 

Since U(-,y) is harmonic, we have 

d i+k U(x,y) dWf a , y (\x-y\ 2 + e ) 



dz a 



z=y 



d e P(x,z) / d k T(y-z) 
1 3D dxP \ dy a 



Idi 



fa, y (\y- z \ +£o))da(z) 



J G(x, z)A z f a>y (\y - z\ 2 + eo) dz 



dx 13 

= Ji(x,y) + J 2 (x,y) 

Notice that the term involving integration over the interior comes from Green's 
identity — since the integrand has a singularity. 

Since dD is C k °' e we have (from the maximum principle) that 



d^U(x,z) 




d^U(z,x) 


dz~< 




dz~< 



< C 7 , e ||r(a; - OllcN(a£») < C na 6(x) 



for all |7| < k . Moreover, since again U(x, •) is harmonic, we have 



-\j\-n+2 



(2.11) 



dW + WG(x,. 



dx^dz-y 



for all I7I < ko and any (5. 

Since P(-,z) (for z £ dD) is harmonic, we have 



dx? 



P(x,z) 



< C e S(x) 



-n-i+l 



for all z e dD. Thus 



r 6(rnY +k+n 6(v) n+2k 



< a 



6{y) 2k+n 



dD 6(y) 



n+k- 



(|^ — y| 2 + e ) 
-da(z) 



< Ck,e- 



BOUNDARY VALUE PROBLEMS FOR THE INHOMOGENEOUS LAPLACE EQUATION 

Also, for R = \x - y |/2 > fixed, 



7 



\M x o>y)\ ^ 



if 



\ G(x,z)A z f ay (\y - z\ 2 + e )dz 

JD\B(v,R) 



+ 



dx 13 JD\B(y,R) 

&G{x Q ,z) 



X=Xq 



I 

Jd 



DnB(y,R) Ox 13 

J2i(x ,y) + J 2 2(x ,y). 



Azfa,y(\y ~ A 2 + e o) dz 



Now 



Jn{x Q ,y) < 



T(x - z)AJ ay (\y - z\ 2 + e ) dz 

OXP JD\B(y,R) 



X=XQ 



+ / e 5(x )- n - e+2 (\z -x \ + R)- n ~ k dz 

JD\B(y,R) 

< J2ii + C k Ax -y\- n - k - e+2 , 



where J 2 n is defined by the last inequality. We see that 



J2ii(x ,y) = 



< 



Ox 13 JD\B(y,R) 



X = XQ 



+ 



+ 



dxP Jd(D\B(y,R)) 

d 



T(x - z)A z f a; y(\y - z\ 2 + eo) dz 

F(x - z)D u(z) f a , y (\y - z\ 2 + eo) da(z) 



X = XQ 



C l r 
dx 13 Jd( 



dxt 3 Jd(D\B(y,R)) 

d e f a , y {x,y) 



D u{z) T(x - z)f aty (\y - z\ 2 + eo) do(z)\ x -_ 



-X{) 



dx? 

< C k/ \x -y\ 



X=XQ 

-n-k-e+2 



by standard arguments (since k < k , dD is C ko,e with k > 2, e > and by the fact 
that \v(z) ■ (x - z)\ < C(\x - z\ 2 + 5(x) for all z G dD and x G D). [Note that the 
last term in the penultimate line comes from Green's theorem.] 
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Now we consider J2 2 (xo, y)- Let Q = D D B(y, R). After applying the divergence 
theorem many times, we have 



< 



d e G(x ,z) 

d e G(x ,z) d k 
< ! O.r 1 dz a 
d e+k G(x , z) 



A z f a; y(\y - z\ 2 + e )dz 



A z f(\z-y\ 2 + e )dz 



+ 



AJ(\z-y\ 2 + e )dz 
d e+ ^G(x ,z) 



dz a -t 



In dx?dz a 

L s i w 

jm h\<k 

= J22i(x ,y) + J 2 22(x ,y). 
Now, since A z f(\y — z\ 2 ) = for all z ^ y, we have 



U -AJ(\z-y\ 2 + e )d ( r(z) 



(2.12). A z f(\y-z\ 2 + e ] 



n(n - 2) e 5(rr ) n+fc+£ %) 2fc+n 



Moreover, we have 



£>^W(|y-^ + e ) 



(| z _ y |2 + eo) n/2 | z _ y |n 

<5(x ) fe+ ^ n %) n+2fe 



5; C n . 



(\z - y\ + e ) n +h\ 



Applying (2.12), we find that 

^222(^0, y) 



< Ct, k Jj(x y 



-fc-n+2 



5{x o y +n+k 5{y) 2k+n 



n+k 



da(z) 



Finally, we estimate ^21(^0,2/)- Observe that 

J22i(x ,y) 



d k+e G(x , z) 



< J C k , n 5(x )' 



A z f(\y-z\ 2 + e )dz 



5{x ) k+n+e 5{y) n+2k 



■k-i-n+2 



i\y-A+^y 



dz 



< C n , k eS +k+e log 



< C n k- 



e 0, 
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Combining the above estimates, we have 

d k+e U(x ,y) d e f a (\x -y\ 2 + e ) 



dx? 



< C n , k (S(y) 



\x 



y\) 



-n-k-i+2 



Q x PQ y a 

Using this estimate and the definition of f a ,y(xo), we conclude that the proof of 
Theorem 2.2 is complete. D 

Corollary 2.3. Let A > 2 be a non-integer and let D be a bounded domain with A 
boundary. Then for all < k < [A] and any cube Q in M. n with 2Q C D, we have 



(2.13) 



d 2 U(x, y) 



dxjdxj 



v d a d 2 U(x,x ) 



(y-xo) a 



< a 



p,n,\" 



X 



- x \ n+k W 



for allx E D\2Q and y,x eQ and \Q\ = 5 n . Where fc(A) = k + 1 if k < [A], and 
fc(A) = A if k = [A] 

Proof. This follows directly from Theorem 2.2 (with a suitable modification for 
fractional derivatives) and from Taylor's theorem. 

3. Reduction of Theorem 1.2 

Let / be a measurable function on D and let A and q be positive numbers. Let 
us recall the following definition for a p-atom. Let a be a bounded function in M, n ; we 
say that a is an h p (R n ) atom if a is supported in a cube Q with (/ R „ \a(x)\ 2 dx) 1 ^ 2 < 
|q|1/2-i/p anc j either (i) \Q\ > 1 or (ii) \Q\ < 1 and, for each a = (a±, ■ ■ ■ ,a n ) with 
\a\ < [n(l/p— 1)] , we have 

x a a(x)dx = 0. 



Jq 



where all ctj > and \a\ — Yh=i a «- 

In order to prove Theorem 1.2, we recall the following theorem that is formulated 
from the statement of Theorem 3.2 in [CKS] and its proof. 

THEOREM 3.1. Let D C W 1 be a bounded Lipschitz domain and let < p < I. 
Then f G h p z (D) if and only if f has an atomic decomposition, with h p (M. n ) atoms 
whose supports lie in D, i.e., 

f = 2 X Q a Q 

QcD 

with 2Q fl dD = if the diameter of Q is small (i.e. < 1), and 

E |AqI p <oo. 

QcD 

By Theorem 3.1, it is easy to verify that Theorem 1.2 holds if and only if it holds 
for / an /i p -atom. More precisely, the proof of Theorem 1.2 can be reduced to proving 
the following theorem: 
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THEOREM 3.2. Let < p < I. Let D be a bounded domain in R™ with C a{p ' e) 
boundary. Then 



d 2 G{o) 



dxjdxi 



+ 



h v r (D) 



d 2 N{o) 



dxjdxi 



h v r {D) 



< C p , n>e (D). 



holds for all p-atoms a with support Q satisfying 2Q C D. 

We shall separate the proof of Theorem 3.2 and Theorem 1.3 into two cases, 
which are given in Sections 4 and 5 respectively. 



4. The Dirichlet Problem 



In this section, we shall prove Theorem 3.2 with the Dirichlet boundary condition. 
In other words, we shall prove Theorem 4.1. 

THEOREM 4.1. Let < p < 1 and let D be a bounded domain in R™ with C a{p ' e) 
boundary. Then 

d 2 G(a) 



(4.1) 



dxjdxi 



h p r (D) 



< C p , n , e (D). 



holds for all p-atoms with support Q satisfying 3Q C D and for any e > 0. 
Since 

G(a)(x) = j D (r(x -y)- U(y, x))a(y) dy, 

it is clear that 

d 2 



dxidxj Jd 



T(y - x)a(y) dy 



< C 



hl(D) 



from the results in [ST1] for Hardy spaces in R n . Let a be an /i p -atom. Without loss 
of generality, we may assume that G D and 

supp(a) C {y e D : \y\ < 5} = B(0, 5) C D. 

If the atom a is supported in a cube Q C D with diameter greater than 1, then 
a E L 2 (R n ). Then we have that (4.1) holds by using a modified version of the 
following result of Jerison and Kenig in [JeK] and Gilbarg and Trudinger [GT]. 

THEOREM 4.2. Let D be a bounded domain in R n with C 2 boundary. If u is a 
solution of (1.1), then 

\\ U \\w% +ct {D) < C p , a \\f\\ W P( D ). 

for all 1 < p < oo and — 1 < a < 0. 
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Note, in passing, that there is certainly an analogous version of Theorem 4.2 for the 
Neumann problem. 

Now we assume that a is a classical atom with support in a cube Q with small 
diameter with 2Q n dD = 0. Let the center of Q be x = 0. For each x G D, we let 



H(a)(x) = ^—J D U(x,y)a(y)dy. 



dxidxj 



Since U(x,y) is harmonic in both x and y, we have that H(x) is harmonic in x. We 
choose G Cq°(B(0, 1)) to be a non-negative radial function such that J R „ <j)(x)dx = 
1. Let x E D and (i(x) = dis(x, -D c ). Since iJ(a) is harmonic and is radial, the 
mean value property of harmonic function shows that 



H(a) + (x) = sup 

0<t<d{x) 



D 



H(a)(y)(j) t (x - y)dy 



\H(a)(x)\, 



where (j> t — t n (j)(t 1 x). Now 



H(a) + (xfdx = / \H(a)(x)\ p dx + / \H(a)(x)\ p dx = Ua) + I 2 (a) 

l) JB(0,4<5) JD\B(0,4S) 



By (1.3), we have 



/ \ P/2 / \ (2-p)/2 

A (a) < f / |#(a)(x)| 2 cfe ) I / dx) 



/ \JB(0,4S) J 
\ P/2 



< C7 p [ / |a(x)| 2 dx |B(0, 45) | (2 - p)/2 

'5(0,4(5) / 

< C7 p |5(0,45)| ( - 2/p+1)p/2 |fi(0,45)| (2 - p)/2 
= C7 p . 



Next we estimate /2(a). By Corollary 2.3, we have 



9^9^ |c ^ p <9t/ a <9x;<9xj 



< C n 



£n(l/p-l)+e/2 



\ x - Xo |n+n(l/p-l)+e/2 



12 



STEVEN G. KRANTZ AND SONG-YING LI 



for all x e D \ AB(x ,5) and y e B(x ,25) C D, where n p = [n(l/p - 1)]. Thus, 
since the center of Q is xo = 0, for each x <E D \ 73(0, 45) we have 

d 2 U(x, y) 



\H(a)(x) 



q dxidxj 
d 2 U(x ) y) 



a(y) dy 



q \ dxidx 



_ d a d 2 U(x,0) a 



< a 



p.c 



3 |q|<"-p 



\X — X 



n+n(l/p-l)+e/2 
£n(l/p-l)+e/2 



a{y)\dy 



x |n+Ti(l/p-l)+e/2 



— C'r, 



X 



n/p+e/2 ' 



Therefore, if D C B(0, d ), then 



D\B(0,4<5) 



|ff(a)(x)|'(fa < C£ e / 

/■do 

< c v I c n 

P J AS 



6 pe/2 

D\B(0,45) \x\ n+pe / 2 
d S pe/2 



dx 



CI 



fl+pe/2 

2 5 pt ' 2 



dt 



V (45) 
= C(p,n,d ). 
The proof of Theorem 4.1 is thus complete. D 

THEOREM 4.3. Let < p < oo and let < n(l/p— 1) be a non-negative integer. 
Let D be a bounded domain in W 1 with (7 2 +n(i/p-i)+e b ounc [ ar y f or some e > 0. Then 

d 2 G(f) 



(4.2) 



dxidxj 



<C k 



h k r ' p (D) 



k,p,e\\J\\hz' p (D)- 



for any e > 0. 
Proof. Observe that 

<?(/)(*)= / T(x-y)f(y)dy- [ U(x,y)f(y)dy = T(f)(x) + U(f)(x) 
Since / has compact support, it is obvious that 



(4.3) 



v2r (/)| W) ^^ll/ll fc ;^)- 



Next we prove that 



(4.4) 



V 2 W) 
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Now V U(a)(i) is harmonic for any non-negative integer £, so it suffices to show that 
||V^C/(a)||x, P (D) < C p for all < £ < k + 2. We shall prove the case £ = k + 2; the other 
cases are even easier. In order to do this, we need the following Sobolev embedding 
theorem which is a special case of Theorem 2 in [BB] and can also be deduced from 
Theorem 2 in [HPW]: 



and the embedding is continuous. By Theorem 3.1 and the fact that V e U(f) is 
harmonic, it suffices to prove that 



for any pn/(n — pk)-aiom with support Q and 3Q fl dD = (since k < n(l/p — 1), 
hence pn/ (n — pk) < 1.) 
Notice that 



(4.5) 



h k /{D) C hl n/in - pk \D) 



d k+2 U(a) 



dx a 



Lp(D) 




Let x be the center of Q. Then 




We set 



£ p = [n(n — pk)/np — n)] = [n/p — k — n], rj(k) — 1 if k > 0, 77(0) = e. 



Then 



d k+2 U{x,y) 
dx a 



E 



ga d k+2 U(x,X Q ) 



(y-x ) a <C P 



\a\<£ p 



x _ ;ro |(n/p-A:)+»;(fc)+fc 
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for all x E D \ 2Q and y E Q as above. Thus, for any x E D \ 2Q we have 
d k+2 U(a)(x) 



dx c 



d k+2 U(x,y) 



q dx a 
^ d k+2 U(x,y) 



a{y)dy 



dx a 

- pAn ' e J Q \x- 



~ L 7^ — n (y-x f)a(y)dy 



\p\<e p 

fi(n/p-k)-n+rt(k) 



dy? 



x 



k+(n/p—k)+r](k) 



dx a 



\a(y)\dy 



< c pAn M l - {n - pk)/pn - 



fi(n/p—k)—n+ri(k) 
x — x ^k+(n/p-k)+r](k) 



a 



p,k,n,e i _ x \k+(n/p-k)+r)(k) 



for any e > 0. Therefore, if D C B(0, do), then 

\H{a){x)\ p dx 



D\2Q 



< ci 



fipri(k) 



P ' k ' n ' e JD\2Q \X - Xo \pk+n-pk+p V (k) 



dx 



< nP xpv(k)-pv(k) 

— ^p,k,n,e,d ° 

The proof of Theorem 4.3 is therefore complete. D 

5. The Neumann Problem 

In this section, we shall prove Theorems 2.2 and 1.3 with the Neumann boundary 
condition. More precisely, we shall first prove the following theorem. 

THEOREM 5.1. Let < p < 1 and let D be a bounded domain with C a{p ' e) 
boundary. Then 

dN(o) 



(5-1) 



dxidxj 



< C P (D) 



K{D) 



for all h p (R n ) atoms with support Q C 2Q C D. 

Proof. Let a be an h p z atom with support Q and \Q\ = 5 n . If S > 1, then a E L 2 (M"). 
By the version of Theorem 4.2 that holds for the Neumann problem, we have that 
(5.1) holds. 
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Now we assume that a is a classical atom with support Q and 2Q n dD = 0. 
In this case, (4.1) follows from the argument of the proof of Theorem 4.1 and the 
following two results (Theorems 5.2 and 5.3.) 

THEOREM 5.2. Let k > 2 be any positive integer and let D be a bounded domain 
with C k °' e boundary. Then, for any multi-indices a and f3 with \a\ = k < k and 
\P\ — £ < k , we have 

d k ^V(x,y) 
1 ' 1 dx?dy a 
for all x, y G D any e > 0. 



\x-y\ + 5{y)) n + k+i - 2 



Proof. We write f(t) = -t 1 -^ 2 for all t > 0. Then F(x - y) = f(\x - y\ 2 ). Let 
C D = Vol(D)- 1 . Then 



V(x,y) = [ P(x,z)V(z,y)dcr(z) + / G(x,z)C D dz 

JdD JD 



'8Y(x — z) dU(x,z) 
dD \ dv(z) dv(z) 

= I dT( f- Z) V{z,y)da{z) 
JdD dh>{z) 

dU(x ) z) 



\V(z,y)dcr(z) + J G(x,z)C D dz 



dD dv(z) 

- [ dU{ *\ Z) V(z,y)da(z)+ [ G(x,z)C D dz 

JdD V[Z) Jd 

h(x,y) - h(x,y) + h(y) 



It is obvious that h(y) satisfies (5.2) by using estimate (2.8). Now we consider 
I 2 (x, y). Notice that 

JdD Oh>{z) 

+ [ A z U(x,z)V(z,y)dz- I U(x,z)A z V(z,y)dz 
Jd Jd 

= ( U(x,z) dT !: Z ~ y ^ da(z) + 0- ( U(x,z)C D dz 
JdD ou(z) Jd 

f dY(z — y) f 

= / F(x — z) — - . . — dcr(z) — Cd I U(x,z)dz. 
JdD oi>{z) Jd 

Using arguments similar to those in Section 2, we see that I 2 (x,y) satisfies the esti- 
mate (5.2). 

For any xo, y € D, if \x — y\< 45 (y), then it is easy to prove that (5.2) holds by 
replacing (x,y) by (x ,y). Without loss of generality, we may assume that \x — y\ > 
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45(y). We consider the term Ii(x,y). If we let 

U{\z -x\ 2 + e ) = (\x-z\ 2 + 60) 1 ""/ 2 , e = 5{x Q ) n+ " +k 5{y) n+2 \ 



then 



h(x,y) = I dV( f Z) V(z,y)da(z) 
JdD dh>{z) 



We have that 



JdD Oh>{z) 

f ( dY{x-z) df xfi (\z-x\ 2 + e ) \^ ( 
hi(x,y) + In(x,y) 



hi(x,y) = [ f x ,o(\z -x\ 2 + eo) ^-rv- da{z) 

JdD OV\Z) 

+ J AJ xfi (\z - x\ 2 + € )V(z, y) dz 



- / fx,o(\ x ~ z \ 2 + e o) C D dz. 
Jd 



By the definition of f x ^ and by (2.12), one can easily see that In(x ,y) satisfies the 
estimate (5.2). 

Finally, for convenience, we assume that n/2 is a positive integer (otherwise we 
use (n+ l)/2 instead). Then 



{u(z),z-x) (u(z),z-x) \ 

^X^zf (| X _ Z |2 + eo )n/ 2 J^y)^) 



= <" - 2 ) "I L ( \x - zn~ 'A e )^ V) d ^ 

By the definition of eo, one can easily see that Iuixo, y) satisfies the desired estimate 
(5.2). Therefore, combining the above estimates, the proof of Theorem 5.2 is com- 
plete. D 



As a corollary, we have the following result. 
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THEOREM 5.3. Let < p < 1 and let D he a hounded domain with C a ^ 
boundary. Then 



(5.3) 



d 2 V(x,y) 



dxjdxi 



|a|<np 



dy a dxidxj 



< a 



£n(l/p-l)+ 7 (e) 



p,n,e 



X 



— X \ n+n( - 1 / p ~ 1 ^^ 



for all x G D \ 2Q and y, Xo G Q C 2Q C -D and any e > 0, where 7(e) = n(l/p — 
l)-[n(l/p-l)] + e/2<l. 

In conclusion, combining Theorems 4.1 and 5.1, the proof of Theorem 1.2 is 
complete. D 

With the same argument, the Neumann problem in Theorem 1.3 can be proved 
by using Theorem 5.3 and the argument for proving Theorem 4.3. We leave the 
details for the interested reader. The proof of Theorem 1.3 is complete. D 



6. Counterexamples 

In this section, we shall give some examples to show that the hypothesis on 
smoothness of dD in Theorem 1.2 is essentially sharp. We first prove the following 
lemma. 

Lemma 6.1. Let D be a bounded domain in M n with C 1 boundary. If (integration 
against the kernel) q x .q x . is bounded from h p z (D) to K p r (D) then, for any xq near dD 
with r(xo) = S(xo)/C(n), we have 

d 2+n p Q 



(6.1) r(x ) prip ~ n(1 - p) 



D\B(x ,r{x )) 



dxidxjdy 



(x,x ) 



dx < C p 



for all multi-indices a with \a\ = n p = [n(l/p — 1)] + 1 and p < 1/2. 

Proof. Let Xq G D be near the boundary and let t(xq) = 5(x )/C(n) with 
C(n) » 1 a constant to be chosen so that there is cube Q(x ) with center at 
x with B(x , r(x )) C 2Q(x ) C D. We consider the function 

(6.2) gM=Mx)= ^^^y 

[Here is a radial bump function as usual.] For each < p « 1, we define 

(6.3) a{x)=r{x Q T^l^ d -^^. 
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It is easy to show that a is a p-atom with support in B(xo,r(xo)) C 2Q C D. Thus, 
for all x G D \ B(x ,r(x )), since G(x, •) is harmonic in B(xo,r(xo)), we have 

d 2 G(x,y) 
d dxj dxj 

d 2+n ?G(x,y) 



±r(x u )"°-' 1/ "-» ffl+ "'' G(3: ' i/) 



d dxidxjdy 
2+n "G(x, 
dxidxjdy 



y=xo 



for all multi-indices a with |a| = n p . Therefore the fact d d x .Q x . maps h p z (D) to h p r (D) 

boundedly implies that (6.1) holds, and the proof of Lemma 6.1 is complete. D 

Let B 2 be the unit disc in M. 2 and let ip be a conformal map from B 2 to some 
domain DCC« R 2 .. We write D = ip(B2), and ip(z) = ip~ 1 (z). Then the Green's 
function for —A on D is 

G D {z,w) = — log- 



2tt \<p(z) - <p(w)\ 
for all z, w E D. Now we have 

Proposition 6.2. With notation above, 



dw 2 pdz 2 dw 2 p \(p(z) — (p(w) J dw 2 p \((p(z) — (p(w)y 

for all z^w, where 2 P = [2(1/ p - 1)] + I. 

Proof. We calculate that 

d 2+2 pG(z,w) 



4:71- 



dw 2 pdz 2 

Q2+2 P 

dz 2 dw 2 p ( lQ S 1 1 - ^Mw)\ 2 - log \(p( Z ) - <f(w)\ 2 ) 



d 2 p ( -<p"(zMw) <p'(z)Mw) 



dw 2 p VI — (f(z)ip(w) (1 — ip(z)ip(w)) 
tp»{z) + V '(z) 2 \ 



2 



2 



<p(z)-<p(w) (<p(z) - <p(w)) 
,n , 9 2 p f 1 \ , 2 d 2 » 



(9w 2 p VV 9 !- 2 ) — V 9 ^) / 9-U7 2 p V(^(-2) — (f(w)) 2 J 

and the proof is complete. D 

Without loss of generality, we may assume henceforth that 2(1/ p — 1) > 2 is an 
integer. 
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Proposition 6.3. Suppose that if) G C /3 (B 2 ) (the same for tp) with f3 = 2(l/p — 1) — e 
(for this last, it is sufficient that dD he C ,2 ( 1,/p_1 ) _e/ ' 2 . Assume that g ® q x induces a 
bounded operator from h p z (D) to h%(D). Then 

r(x y\^\x )\ p C m) < C p (l + |M| A/3 ) 2 ", 

where 

2<p'(z) 2 <p"(z) p 



C<p,p,x = 

Proof. Now 

47T 



D\B(x ,r(x )) 

d 2+2 ?G(z,w] 



(<p(z) - (f(x )) 3 ((p(z) - <p(x )) 



dA(z). 



> 



dz 2 dw 2 p 

2ip'(z) 2 ip^\w) <p"(z)<p^(w) 



(ip(z) - ip(w)) 3 (ip(z) - (f(w)y 



_ 2 g 2 C p (l + |^| c2p - 2(5) ) 2 - C p (l + y\\ Ag(B) )5(w)- 



> 



\(p(z) - <p(w)\' 



2 V '(z) 2 



((p(z) - <p(w)) 



\<p(z) - <p(w)\* 
<p"(z) 



\<p(z) - ip(w)\ 2 + 2 ? 5(wY\ip(z) - y?H| 4 ' 
We conclude that 



y? (2p) (u>) 



\ip(z) - ^{w)\< 



2(p'(z 



(ip(z) -<p(w)) 



tp"(z) 



< 4vr 



d 2+2 "G(z,w) 



dz 2 dw 2 p 



C r 



p\m\Ap 



\<p(z) — (p(w)\ 2+2 p 5(w) e \ip(z) — (p(w)\' 



Also 



and 



ix yJ L 



1 C p r{x ) p 



ix ° r J L 



D\B(x ,r(x )) \ip(z) - ^(x )\P( 2+n p) 
1 



dA(z) < 



C n 



r(x ) 2 +P- 2 
dA(z) < C p r(x )P-P € < C p . 



(t) 



lD\B(x ,r(x )) 5{Z)P^{Z) - ^{X )\ 4 P 

Since 2(l/p — 1) is integer, we have 

P 2 P - 2(1 -p)= p(2(l/p - 1) + 1) - 2(1 -p)=p. 

Combining this and all those estimates with Lemma 6.1 and Propositions 6.2, the 
proof of the proposition is complete. D 
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Proposition 6.4. For any e > 0, there is a bounded domain D in M? with C 2 ^ 1 ^ ^ e 
boundary and the operator induced by dx .g x . is not bounded from h p z (D) to h p (D). 

Proof. Let 

^(z) = z + v (l - zjWHH 

with < i] = r](p) « 1 sufficiently small, depending on p. Clearly, 

il>\z) = 1 + n(2(l/p - 1) - e)(l - zf' p -^ 

It is easy to see that ip is a conformal map from B 2 onto ip(B 2 ) provided that < rj = 
n(p) is sufficiently small and ip G Ap(B 2 ) where (3 = 2(1 /p — 1) — e. Let D = ip(B 2 ) 
and <f(w) = ^^(w). It is clear that 



C 



1 . S) > 1/C > 



(p,p,ip(l—8) 

(refer to equation (f)). Now, since 2 P = 2(1/ p — 1) + 1, we have 



(2 P )# 



-e-1 



where c Pit ^ 0. [Note that superscripts in parentheses denote derivatives.] We let 
Ho = 1 — 5. Then 

v 



^ 2p) (y ) 



= M5 p 5 



-ep-p 



= C 



p.c 



Pfi-pe 



OO 



as 5 — > + , since </? has boundary behavior similar to that of tp. Thus 



( 2 *>(^(l-<f)) 



Seeking a contradiction, we suppose that 



d 2 G 



induces a bounded operator from 



h p z (D) to h^D). Then, by Propositions 6.2 and 6.3, we have 



for any < 5 << 1. This is a contradiction as 5 — > + . Thus the operator induced 



by 



a 2 G 



is not bounded from h p (D) to h p (D), and the proof is complete. D 
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